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Introduction 


¢ Many engine power-loss events reported since the 1990's 

¢ Ice crystals entering the engine core — Mason et al. 

¢ Ingestion of ice into engine studied at NASA PSL and elsewhere 

¢ Observed environmental conditions changed with cloud activation 
— Gas temperature change 
— Humidity change 

¢ Hypothesis: Thermal interaction between air and cloud 

¢ Building on previously written model to simulate PSL 

¢ Objective: Understand the air - cloud interactions in PSL tunnel 
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General Description of Model 


¢ Model Simulates PSL icing tunnel 


S 


— Air and cloud conservation equations (mass, energy) fully coupled 


— A\tr is treated as ideal compressible gas 


— lsentropic equations used to solve pi, Vairs Tair P 
— Air and particle flow are steady and one dimensional 
— Temperature is uniform within the perfectly spherical particle 


— Full particle size distributions used 
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PSL Geometry and Capabilities 


Tunnel Capability 
¢ Freeze out liquid cloud 


¢ 12 parameters can be varied 


— PV, Tair Twates RH, MVD, TWC, 
Water Type, Nozzle Pattern... 
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Tunnel Controllability 
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Supersaturation and Condensation 
¢ Vapor saturation can be exceed for certain conditions 
¢ Condense on cloud particles through diffusion not sufficient 
¢ Supersaturated? Condense? Combination? 


¢ 2 type of condensation 
— Homogeneous - RH >> 100% (very clean air) 
— Heterogeneous - RH >100% (nucleation / seeding) 


¢ Nature ~ 101% RH 


RH, = 127% 
U = 135 m/s 


RH, = 76% 
U=3 m/s 
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Condensation Cloud Experiments 
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Aerosol Particulates Background 


¢ Organic and inorganic in composition 
¢ Size distribution from 0.003 um to 2.5 um 


¢ # density variations 
— 3,100/cms (Alps) 
— 100,00/cms (city background) 
— Diurnal variation (peak traffic hours) 
— Seasonal variation (heating in winter) 


¢ Aerosol particulates considered in condensation 
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Aerosol Condensation Subroutine 


¢ Implemented only when RH>100% 
¢ Treat aerosol like any other water droplet / ice particle 


¢ Initial # Density: 22,000/cms (Pittsburg, PA paper) 
¢ Initial Size: 0.04 um (Pittsburg, PA paper) 

¢ Initial Velocity: 99.99% of air velocity 

¢ Initial Temperature: Two 

¢ Twb >0 °C : Condense as liquid 

¢ Twb <=0 °C : Deposit as ice 


¢ Effects of charged particles neglected 
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Model Formulation - Algorithm 


¢ Writtenin MATLAB version R2016b 
¢ Solves conservation differential equations using built-in ODE45 solver 


¢ Numerical relative and absolute convergence tolerance of 10° 
¢ Mass transferred between the gas and particle(s) balanced to 10°'® 


¢ Energy transferred between the gas and particle(s) balanced to 104 
— Physicalaccuracy dependent on accuracy of property values (C,, 
Lieat etc.) 
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Supersaturation Simulation Profiles 
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Supersaturation Simulation Comparisons 


(w= mass mixing ratio) 


pai2ial;a{s}o6 {7 {eo | wo | a 
Cond | Spray | 79; | Tecate} RHoi |RHs,¢,calce}] “100%RH j , We sim,none| Ve,sim,aero 
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Experiment Configurations 
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Tomography — Icing Cloud Spread 
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Experimental Test Conditions for 4 RH Sweeps 


¢ Varied Parameters 
— RH ;=0% to 60% 
— MVD;= 15 um or 50 um 
— U,= 85 m/s and 135 m/s 


¢ Constant Parameters 
— Ty; =7.2 °C 
— Po; = 44.6 kPa 
— TWC;=7.0 g/m? 


¢ Twb Ranges 
— Twbp;= -6.9°C (0% RH) 
— Twhoj; = +2.4 °C (60% RH) 
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Plenum RH Sweeps - AHumidity 
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Plenum RH Sweeps - AT,,, 
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Wet-Bulb Temperature Change, °C 


Wet-Bulb Temperature Change, °C 
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Plenum RH Sweeps — Melt Fraction 
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Plenum RH Sweeps - TWC 
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Plenum RH Sweeps - MVD 
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Summary 


¢ Model written to understand Air - Cloud interactions in PSL 
¢ Aerosol Condensation implemented for better accuracy 
¢ Model over-predicts amount of evaporation (AT,;,, A Hum) 


— Correcttrend for varying RH 
¢ Smaller Two changes, important to determine cloud phase 
¢ Good agreement for melt ratio 
¢ TWC and MVD comparisons suggest 2D effects 
¢ 1D model will not capture 2D cloud movement 
¢ Provides useful predictions even as 1D 


— Model guided development of test matrix for fundamental ICI tests 
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Simulation Results — Aerosol Parametric Analysis 


Test Conditions 
Ty; = 10.0 0C U, = 135 m/s 
Po; = 78.2 kPa MVD, = 15 um 


RH); =77% TWC, =7.1 9/m3 
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